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Abstract 

We present Perl Algorithm to Compute continuum 
and Equivalent Widths (pacce). Wc describe the 
methods used in the computations and the require- 
ments for its usage. We compare the measurements 
made with PACCE and "manual" ones made using 
IRAF splot task. These tests show that for SSP mod- 
els the equivalent widths strengths are very similar 
(differences <0.2A) for both measurements. In real 
stellar spectra, the correlation between both values is 
still very good, but with differences of up to 0.5 A. 
PACCE is also able to determine mean continuum 
and continuum at line center values, which are help- 
ful in stellar population studies. In addition, it is 
also able to compute the uncertainties in the equiv- 
alent widths using photon statistics. The code is 
made available for the community through the web at 
[http://www.if.ufrgs.br/~riffel/software.html . 

Keywords Methods: data analysis; Line: Equivalent 
Widths; Techniques: Spectroscopy 



1 Introduction 

The equivalent widths {W\) of absorption lines, ob- 
served in the spectrum of astronomical sources can be 
seen as compressed, but highly informative, represen- 
tation of the whole spectrum. For example, the W\ 
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of the absorption lines observed in galaxies spectra re- 
veals insights about their stellar populations, like the 
ages and metallicities of the stars which dominates the 



light oi the host galaxy 
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ing the spectrum of a star, it is possible to make use 
of the W\ of absorption lines to determine directly 
the fundamental atmospheric parameters such as: sur- 
face gravity (logg), effective temperature (Teff) and the 
chemical abundances of many elements (see for exampl e 
Gonzalez fc Lambertl[l996i iFeltzing fc GonzalezlboOll) . 



However, the price to be paid when using the pow- 
erful informations contained in the W\ is the long time 
needed to make a reliable measurement of this observ- 
ables. Commonly, the W\ are measured using inter- 
active routines like splot provided by t he IRAILj tea m 



or wi th independent codes like LINER (jPogge fc Owen 
19931) . Both softwares are "hand operated" , which 
means the user need to look for the line limits and con- 
tinuum points in the spectrum. The next step, is to 
mark them "manually" . This procedure is very time- 
consuming and introduce many uncertainties which are 
propagated to a posterior analyses of the quantities in- 
volving the Wx measurements. 

With the growing of spectral surveys (e.g. Sloan 
Digital Sky Survey), it becomes necessary to acceler- 
ate and automate some process such as the analysis of 
stellar populations of galaxies, as well as the determi- 
nation of fundamental atmospheric parameters of indi- 
vidual stars. In order to help in such task we present in 
this paper a new automatiqj code: PACCE: Perl Algo- 



^IRAF is distributed by National Optical Astronomy Observato- 
ries, operated by the Association of Universities for Research in 
Astronomy, Inc., under contract with the National Science Foun- 
dation, U.S.A. 

^Also interactive, as the input parameters are easy to be changed. 
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rithm to Compute Continuum and Equivalent Widths. 
This software, written in PERL, can be used to com- 
pute the W\ of absorption Hnes as well as to determine 
continuum points, being very helpful to perform stellar 
population sy nthes i s follo wing , for example the met hod 
developed bv lBical (Il988l ) and lschmitt et all |l996l) . 

This paper is structured as follows: In Sec. [5] we 
describe the system requirements as well as the numer- 
ical procedures behind the code. The input parameters 
and the outputs of the code are discussed in Sec. |3l A 
comparison with the measures made whit pacce and 
"hand-made" measurements are presented in Sec. SI 
The final remarks are made in Sec. [3 



2 The code 

The idea behind PACCE is to reproduce the "manual" 
procedure used to measure the W\ of absorption lines 
in a spectrum, as well as to measure mean continuum 
fluxes in defined regions and compute the continuum 
value at line center. In addition, using the same in- 
puts, it does exactly reproduce the measured values 
being "user" independent (e.g. the uncertainties intro- 
duced by the user in "hand operated" procedures are 
removed), and thus allowing for a better comparison 
between Wx measured by different users. 

PACCE was written in Perl, allowing anyone to 
use it without having any problems with software li- 
censes. It is freely distributed under GNU General 
Public Licenstll(GLP). All the libraries used in the 
code are also free and distributed under GLP license. 
PACCE source code can be freely downloaded from 

All re- 



http: / /www. if. ufrgs.br/~riffel/software. html 



quirements to run PACCE can easily be installed in any 
linux machine (e.g. apt-get, synaptic or yum), they are: 

• Perl 



(http://www.perl.org); 



Perl's "Math::Derivative" package 
( |http ://search. cpan.org); 
Perl's "Math:: Spline" package 



(http://search.cpan.org 



• Gnuplot 

(www.gnuplot.info) . 

In addition, we call the attention to the fact that Perl 
is installed as default in any linux flavor and PACCE can 
easily be converted to run under Microsoft Windows 
system (without plots). 
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Fig. 1 Fluxogram describing the use of PACCE . 

2.1 Calculations: Equivalent Widths and Mean 
Continuum Fluxes 

2.1.1 Equivalent Widths 

In general, absorption feature indices are composed by 
measurements of relative flux in a central wavelength 
interval corresponding to the absorption feature con- 
sidered (line limits A; and A.^ in Fig. [2]) and two con- 
tinuum sidebands passband regions (spectrum ranges - 
dots - in the boxes of Fig. [2]). Such sidebands provide 
a reference level (pseudo-continuum, solid line Fig. [2]) 
from which t he strength of the a, bsorption feature is 
evaluated (see IWorthev et al. 1994 for details) . PACCE 
computes the pseudo-continuum, Fc{X), in three ways: 

(i) a linear regression is computed using all the contin- 
uum points in the passband regions and a straight line 
(y=ax+b) is computed using the regression coefficients; 

(ii) a straight line is drawn connecting the mid-points 
of the flanking passband continuum regions; (iii) as a 
cubic splin^. The form in which the pseudo-continuum 
is adjusted is chosen by the user (see Sec. [3]). 

Considering Fi{X) as observed flux per unit wave- 
length an Wx is then: 



Wx 



1 - 



dX 



(1) 



^http:/ /www. gnu.org/licenses/gpl. html 



■^Such option is useful when measuring Wx of absorption lines 
located close to emission lines, and thus, only points - not regions- 
of continuum free from emission / absorption arc available 
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Fig. 2 Sketch of the definition of a hypothetical index. 
Sidebands continuum regions are represented by boxes. 
Solid line is the pseudo-continuum. Vertical lines represent 
the line limits A; and \u- Shaded region (A2) is the area 
below the line limits and filled (Ai) region is the area of the 
line. Dotted line (red) represents the spectrum. 



Considering Fig. [2] one can write the area A2 as: 
A2= r F,{\)d\ (2) 

and, thus, the Ai is the absorbed flux between A; and 
A„. Similarly, the area below the pseudo-continuum 
between A; and A^ is; 



Ai+A2=G= I F,{X)dX, 



(3) 



thus, the W\ of the absorption line between A; and A^ 
can be written as: 



(4) 



for mo re details see for example lVollmann. fc Eversberg 
(|2006h . 

In PACCE we compute both areas, A2 and C, using 
the trapezium method. In addition, PACCE does com- 
pute the uncertainties in the equivalent widths, a{W\), 
considering the photon noise statistics pixel by pixel. 
For t his purpose we follow Vollmann. fc EversbersJ 
( 20061 ). assuming that the ratio between C and A2 are 
similar to the norm alized ratio between these qu anti- 
ties (for details see Vollmann. fc Eversberg 20061 ). In 



the case of linear pseudo-continuum adjustment we es- 
timate the signal-to- noise ratio (S/N) as being the ratio 
between the square root of the variance and the mean 
flux of the points in the bandpass interval (i.e. boxes in 
Fig [2]) . In the case of a pseudo-continuum defined with 



a cubic spline the S/N is calculate in the same way as 
in the linear adjustment, but using points in a region 
free from emission/absorption lines defined by the user 
(see Sec. [3]). 

2.1.2 Mean Continuum Fluxes 

Besides the W\ the code also computes the mean con- 
tinuum fluxes {F\, last line of Tab. [1]) and the contin- 
uum at line center. Such measurements are useful, for 
example, to perform ste llar p o pulat ion studies using the 
technique described by Bica ( 1988| ). The ^\ point are 
calculated following the equation: 



^ i=N 

^^^nT. ^a, 



(5) 



where is the flux of each A in the defined interval 
and N is the number of points considered. The errors 
are estimated as being the square root of the variance. 
The continuum at the line center is taken as being the 
pseudo-continuum flux where A 



Ai+A„ 
2 



3 Running PACCE 

To run PACCE the user needs a ASCII, one dimensional, 
spectrum and a ASCII input table containing the line 
definitions. The ASCII spectrum is easily created using 
the wspectext IRAF task. An example of a input table 
is shown in Tab. [T] Note that the continuum points 
used to determine the line pseudo-continuum can be 
defined into three different ways (linear, mid-point or 
spline). The fact that such input table is very easy 
to be edited/created, makes PACCE also an interactive 
tool, allowing for fast changes and tests in W\ and F\ 
measurements . 

Besides the calculations, PACCE outputs some aux- 
iliary files (Fig. [1]). These files are the pseudo- 
continuum data points, the plots showing the regions 
used in the W\ calculations (see Fig. [3]) , as well as the 
Gnuplot commands, which are stored for future use. 



4 Testing PACCE 

In order to test PACCE we perform "ha nd-made" mea- 
surements in a set of iMarastorJ (j2005() simple stellar 
population (SSP) models as well as in a library of ob- 
served stellar optical spectra (jSilva fc Cornelllll992l ). 

In Fig. |3]we show a comparison ^^3600 W\ measured 
using PACCE and ' 'hand-made" measurements with 
IRAF splot task. We also plot the differences between 
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##################################################################################### 

# # 

# Index definitions table example. # 

# # 
##################################################################################### 

# ID Line Limits Left side Cont . Right side Cont . 
LinearC (4847 . 875 , 4876 . 625) [4827 . 875 : 4847 . 875] , [4876 . 625 : 4891 . 625] 
MidPntC (4847 . 875 , 4876 . 625) [4827 . 875 : 4847 . 875] , [4876 . 625 : 4891 . 625] * 

# Note the * in the mid-point continuum adjustment. 

######################### Spline continuum example ################################## 

# ID Line Band Pass Points for a Spline Cont. S/N cont. 
spline (4847 . 875 , 4876 . 625) {4827 . 875 , 4833 . , 4847 . 875 , 4876 . 625 , 4891 . 625}<4880 : 4890> 
##################################################################################### 

# # 

# Intervals for Mean continuum # 

# see Bica (1986) # 
##################################################################################### 
# 

Cont . I 5290 : 5310 , 5303 : 5323 , 5536 : 5556 , 5790 : 5810 , 5812 : 5832 , 5860 : 5880 , 6620 : 6640 I 
Table 1 Example of input table. A # indicates that the line is not used in the computations. Note the * in the mid-point 



continuum adjustment. A full table is available at http://www.if.ufrgs.br/~riffel/software.html 




00 CO CO 



It is even harder to properly measure Wx in ob- 
served spectra than in theoretical ones. To properly 
test our code capability in dealing with real spectra we 
have measured a set of Wx (~ 950W x) in the library 
of obs erved stellar spectra presented by Silva fc Cornelll 
(|l992l) . The results of such test are shown in Fig. [SJ It 
is clear that pacce is able to properly reproduce the 
measurements of splot task. However, a larger disper- 
sion in the differences is observed between both mea- 
surements than when using SSP models. In addition, 
these differences are clearly within the errors and are 
lower than ^O.SA. 



Fig. 3 Example of PACCE output plots using Tab.[T]index 
definitions. Mid-point, linear and spline pseudo-continuum 
adjustments where used in left, center and right plots re- 
spectively. The shown region is around H/3 absorption line. 



the measurements and make a histogram of these dif- 
ferences to give a better idea of the dispersion between 
both sets. It is clear from Fig. Uthat PACCE does re- 
produce very well the splot "manual" Wx values. Note 
that there seems to be a very slight systematic under- 
prcdiction of the Wx computed with PACCE if com- 
pared with those of splot, however, the di fferences are 



<0.05 A. Similar results were obtained by ISousa et al 



1 20071 . see their Fig. 6), thus, suggesting that iraf splot 
measurements may slightly overestimate the Wx values. 



5 Final Remarks 

We present PACCE a Perl algorithm to compute con- 
tinuum and Wx- We describe the method used in the 
computations, as well as the requirements for its use. 

We compare the measurements made with PACCE 
and "manual" ones made using iraf splot task. These 
tests show that for SSP models (i.e. high S/N) the Wx 
values are very similar (differences < 0.2 A). In real stel- 
lar spectra, the correlation between both values is also 
very good, but with differences of up to 0.5 A. However, 
these small differences can be explained by the intrin- 
sic errors in subjectiveness determinations of continuum 
levels caused in "manual" measurements. 

PACCE is also able to determine mean continuum 
and continuum at line center values, which are helpful 



PACCE 
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in stellar population studies. In addition, it is also able 
to compute the a{W\) using photon statistics. 
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Fig. 4 Comparation of pacce and "hand made" iraf SPLOT measurements using iMarastonI ()2005l ) theoretical models. 
Solid line is the equivalence line. At the differences and a histogram of them are shown in bottom and right side, respectively. 
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Fig. 5 Same as Fig. |4]but for observed optical stellar spectra. Error bars are shown. 
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